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The space-borne antimatter experiment PAMELA has recently reported a surprising rise in the 
positron to electron ratio at high energies. It has also recently been found that electromagnetic 
radiative corrections in some cases may boost the gamma-ray yield from supersymmetric dark 
matter annihilations in the galactic halo by up to three or four orders of magnitude, providing 
distinct spectral signatures for indirect dark matter searches to look for. Here, we investigate 
whether the same type of corrections can also lead to sizeable enhancements in the positron yield. 
We find that this is indeed the case, albeit for a smaller region of parameter space than for gamma 
rays; selecting models with a small mass difference between the neutralino and sleptons, like in 
the stau coannihilation region in mSUGRA, the effect becomes more pronounced. The resulting, 
rather hard positron spectrum with a relatively sharp cutoff may potentially fit the rising positron 
ratio measured by the PAMELA satellite. To do so, however, very large "boost factors" have to 
be invoked that are not expected in current models of halo structure. If the predicted cutoff would 
also be confirmed by later PAMELA data or upcoming experiments, one could either assume non- 
thermal production in the early universe or non-standard halo formation to explain such a spectral 
feature as an effect of dark matter annihilation. At the end of the paper, we briefly comment on 
the impact of radiative corrections on other annihilation channels, in particular antiprotons and 
neutrinos. 

PACS numbers: 95.35.-|-d,98.70.Sa,13.40.Ks, ll.30.Pb 



The existence of a sizeable dark matter contribution 
to the total cosmological energy density seems by now 
to be established beyond any reasonable doubt, the most 
recent estimates [l| giving the fraction of cold dark mat- 
ter to the critical density as VLcdm ~ 0.233 ± 0.013. 
Against this background, searches for experimental sig- 
natures that may determine the so far still elusive nature 
of the cosmological dark matter are becoming ever more 
important. 

On the theoretical side, maybe the best motivated, and 
certainly most extensively studied, dark matter candi- 
date is the supersymmetric neutralino (for reviews, see 
0). The methods of detection for this type of parti- 
cle dark matter can be grouped into accelerator searches, 
trying to directly produce dark matter or related new 
particles (the signature of the former usually being miss- 
ing energy), the direct detection of dark matter particles 
scattering off the nuclei of a terrestrial detector, or indi- 
rect detection of particles generated by the annihilation of 
dark matter particles in the Galactic halo or (for neutri- 
nos) in the Sun or Earth. With the LHC soon operating 
and new detectors of liquid noble gases being developed 
for direct detection, aiming to further improve the al- 
ready impressive recent upper limits 0], the near future 
promises very interesting times for the field. 

As far as indirect detection is concerned, the anti- 
matter detection satellite PAMELA has just announced 
its first set of data for cosmic ray antiprotons J| and 
has very recently done so also for positrons [5|. Al- 
though the antiproton data seems to agree with conven- 
tional secondary production by cosmic rays, the positron 



data shows an unexpected rise in the differential ratio 
e+/(e^ -I- e+) above some 7-10 GeV. This interesting 
situation may be further investigated by the PEES bal- 
loon experiment Q and, in particular, the AMS-02 ex- 
periment, if installed on the international space station 
[7| . These experiments could further improve these data, 
both concerning statistics and energy range, and in par- 
ticular investigate whether a return to a "normal" ratio 
at some energy exists - something that is predicted by 
models of dark matter annihilation due to the kinematic 
limit that appears at an energy equal to the dark matter 
particle mass. Of course, it will also be important to rule 
out positron misidentification through proton contami- 
nation in this high-energy range. 

Further information may possibly be obtained from 
the huge IccCUBE Q detector which will soon start to 
look for cosmic neutrinos at the South Pole. For gamma- 
rays, the recently launched GLAST satellite opens up 
a new window to the high-energy universe, for energies 
from below a GeV to about 300 GeV. The sensitivity to 
gamma-rays of even higher energies is, furthermore, ex- 
pected to improve considerably with next generation Air 
Cherenkov Telescopes like the CTA or AGIS USEj. 

Most hkely, in fact, a signal from more than one type 
of experiment will be needed to confirm a dark matter in- 
terpretation of the observed signal and, in the best case, 
to fully identify the particle making up the dark matter; 
it is thus important to realize the complementary nature 
of the methods described above. This is even more true 
since, in all of these cases, the signal searched for may be 
quite weak and dominated by a much larger background. 
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In this context, one should also keep in mind (see, e.g., 
[l^ ) that it may be possible to explain a signal like the 
one recently reported by PAMELA in a more conven- 
tional way, without having to invoke unreasonably strong 
astrophysical sources; for example, a supernova remnant 
of age 10^ years some 100 pc distance from the Sun/Earth 
would both have the appropriate energetics and the right 
energy spectrum to account for the PAMELA results. 

A first assessment of the situation (see, e.g., 
concerning dark matter candidates after the surprising 
PAMELA results seems to indicate that the otherwise 
favoured supersymmetric neutralino cannot explain the 
data. This is because it is a Majorana particle and 
therefore does not give hard positrons directly, due to 
the helicity suppression of light fcrmions in the annihila- 
tion process. The resulting positron spectrum is thus 
expected to be rather soft, in disagreement with the 
PAMELA data, and therefore a Dirac particle, or a spin- 
1 particle like Kaluza-Klcin dark matter IJ] would fit 
better (another proposal put forward in connection with 
the PAMELA data has been minimal dark matter, where 
the combination of a very high dark matter particle mass 
('^10 TeV) and a very efficient enhancement mechanism 
for the annihilation into charged gauge bosons would re- 
sult in the required hard spectrum at low energies [S]). 

However, this simple intuition may prove wrong when 
computing radiative corrections. Gamma-rays, for in- 
stance, have a sharp cutoff 0,113 at an energy equal 
to the dark matter particle's mass, = m^, and, in 
some cases, even prominent line si gna ls from the direct 
annihilation into photons 18, [l9l. l20j. While the ori- 
gin of the first feature is associated with photons di- 
rectly radiated from charged final legs ("final state ra- 
diation" ) , it was recently pointed out that even photons 
radiated from charged virtual particles ("virtual" inter- 
nal bremsstrahlung (IB), or direct emission) can have a 
significant impact on the resulting gamma-ray spectrum, 
leading not only to an even more pronounced cutoff, but 
also to clearly observable bump-like features at slightly 
lower energies [2l|. In fact, these effects generically dom- 
inate the total spectrum at high photon energies, includ- 
ing even the line signals, and may lead to an enhancement 
of the annihilation rate by several orders of magnitude. 
Such a large radiative "correction" can appear since the 
annihilation of neutralinos into lepton pairs is strongly 
helicity suppressed, while for the three-body final state 
containing an additional photon this suppression is cir- 
cumvented 
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With these recent results in mind, the question thus 
naturally arises whether the same effects also have a 
significant impact on, e.g., the yield in positrons - es- 
pecially since the largest enhancement factors appear 
for neutralino annihilations into leptons 



211. 



Let us 

first consider the direct annihilation into positrons. As 
e"'"e~ two-body final states are strongly suppressed, the 
dominant contribution comes always from the process 



XX ~* e"'"e~7, in particular from those diagrams where 
the photon is radiated from a t-channel selectron. Set- 
ting me — > 0, we find for the differential annihilation rate 
into positrons: 
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where X = Ef.+ /m^^ I^r.l 

is the coupling between neutralino, electron and right 
handed (left-handed) selectron. In the corresponding 
limit, this reproduces the result found in for photino 
annihilation. Integrating Eq. ([T]) gives 
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where Li2(2:) = J2T=i •^'^/fe^ is the dilogarithm. The di- 
rect annihilation of neutralinos (with small galactic ve- 
locities v) into positrons is thus only suppressed by a 
factor of (a/vr) and not, as often quoted, by the much 
smaller factor of (mg/m^) that is connected with two- 
body final states. In the above expression, the highest 
annihilation rate is obtained in the limit (ir^l ^ 1 where 
the selectrons are degenerate in mass with the neutralino: 
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Positrons may also be produced in the decay of other 
annihilation products. The number dN^^ / dx of such sec- 
ondary positrons per annihilation into the corresponding 
final state / can be simulated with Monte Carlo event 
generators like PYTHIA [2^. For two-body final states 
XX, we use the tabulated values contained in Dark- 
SUSY [i^] that were obtained through a large number 
of PYTHIA runs. For three-body final states containing 
a photon, the positron yield is approximately given by 
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FIG. 1: The solid line gives the total number of positrons 
per neutralino pair annihilation and positron energy for the 
benchmark model BM3 of [2l[ (m^ — 233 GeV, mg — 
240 GeV). Shown separately is the same quantity without 
radiative corrections (dotted line) and, on top of this, only 
the e^e~7 final states (dashed line). 



where dN^^ /dx is the (two-body final state) positron 
multiplicity dN^j_^ / dx that results from the annihilation 
of two dark matter particles with mass Ex- The ana- 
lytically obtained expressions for dN^^'^ /dEx are too 
lengthy to reproduce here but have been fully imple- 
mented in the current public release of DarkSUSY [2J| 
(for light fermions, of course, the same functional form 
as in Eq. ^ is recovered). When compared to gamma 
rays from the corresponding channel, the positron con- 
tribution ^ to the total spectrum is considerably less 
pronounced at the observationally most relevant energies 
near the cutoff since part of the energy is taken away 
by the photon; the fact that positrons are not the only 
decay products induces a further kinematical suppres- 
sion at high energies. On general grounds, we there- 
fore cannot expect large radiative corrections to the yield 
in secondary positrons - even in situations where large 
gamma-ray contributions are found (as, e.g., for heavy 
neutralino annihilation into W~^W~ [25|). An exception 
to this conclusion could only occur in a situation where 
the annihilation rate into the three-body final state is 
many times lar ger than for the two-body final state. As 
pointed out in [21|. this is indeed possible for lepton fi- 
nal states in the stau-coannihilation region of mSUGRA. 
However, as also the annihilation into e~^e~j is usually 
greatly enhanced in this region, it is, rather, the latter 
contribution that dominates in this case. 

For illustration, we show in Fig. [1] the effect of radia- 
tive corrections on the positron yield for a typical model 
in the mSUGRA coannihilation region (introduced as 
benchmark point BM3 in [IH), which is characterized 
by small mass differences between neutralino and slep- 
tons. A spectacular boost in the positron yield can be 
observed, leading to an extremely pronounced cutoff at 



E^+ = rriy^. As anticipated, this is mainly due to primary 
positrons, following the distribution ([T]), but at smaller 
energies the effect of radiative corrections becomes also 
visible for other decay channels (mainly fi'^fi~). 

The propagation of charged particles is influenced by 
magnetic fields residing in the Milky Way which, in con- 
trast to the case of photons, tend to erase clearly pro- 
nounced spectral features. To be able to compare our 
results with the cosmic ray positron spectrum as mea- 
sured at the top of the atmosphere (TOA) or, in the case 
of the recent PAMELA data, in space, we adopt the stan- 
dard assumption of randomly distributed magnetic fields, 
in which case the determination of the positron flux at a 
given galactic position boils down to solving a diffusion 
equation (for more details on the procedure we follow, 
see d^). We assume an NEW profile 27 1 for the dark 
matter distribution in the galactic halo, but allow for an 
additional "boost factor" to account for the effect of dark 
matter substructure. 

Let us now quantify our general expectations outlined 
above and try to assess the general importance of IB ef- 
fects on the positron yield. Eor that purpose, we consider 
a scan (based on the work in [2^) over the mSUGRA pa- 
rameter space, keeping only models that feature the right 
relic density and pass all current collider bounds. In this 
setup, as mentioned before, the stau coannihilation re- 
gion is characterized by light leptons; since, at the same 
time, the total annihilation cross section (today, i.e. for 
w — > 0) is very small, we expect rather large enhance- 
ments in the positron yield in this case. In the left panel 
of Fig. [21 indeed, we clearly observe the expected strong 
correlation between the rrie-m^ mass splitting and the re- 
sulting enhancement in the positron flux due to radiative 
corrections; outside the coannihilation region, no sizeable 
flux enhancements are encountered. 

In mSUGRA, however, in the presence of light selec- 
trons, also the other sleptons have to be light. Hence, 
to investigate the situation in more general terms, we 
have set up a low-energy phenomenological scan with 9 
free parameters. We start with the usual parameters of 
MSSM-7: the Higgsino mass parameter /x, the gaugino 
mass parameter M2, a common sfermion mass scale toq, 
the ratio of the Higgs vacuum expectation values tan/3, 
the trilinear couplings in the third generation At and Ab, 
and the mass of the CP-odd Higgs boson, tua- In ad- 
dition to these, we add a selectron mass parameter that 
goes into the selectron entries in the mass parameters of 
the soft SUSY-breaking Lagrangian; furthermore, we re- 
lax the GUT condition on Mi and M2 by allowing Afi 
to be varied freely. We then generate models by vary- 
ing these 9 parameters of our MSSM-9 model between 
generous bounds, focusing on models with light selec- 
trons and, again, keeping only models that feature the 
right relic density and pass all current collider bounds. 
In the right panel of Fig. we plot the enhancements 
from IB for these models; the effect of introducing more 
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FIG. 2: Scan over mSUGRA (left) and MSSM-9 (right) models that shows the enhancement in the positron flux (at = 
0.9 mx) due to radiative corrections vs. the mass splitting between the lightest selectron and the neutralino, S = {mg — m^) /m^. 
Also indicated in this figure are the benchmark model BM3 from [2f[ and a further benchmark model BM5' as introduced in 
the text. 



free parameters as compared to mSUGRA is mainly that 
the total annihilation cross section is not anymore closely 
linked to the me-m^ mass splitting by the relic density 
requirement. As a result, the enhancements in the total 
positron flux can be both considerably larger and smaller 
than in the mSUGRA case, depending on the total an- 
nihilation rate to lowest order; the main contribution to 
the flux enhancement, however, is in any case found from 
the e+e^7 channel. 

In Fig. [3l we plot the resulting flux ratio e+/ (e+ + e") 
from neutralino annihilations for both BM3 and a point 
BM5' in the MSSM-9 parameter space (with = 
132 GeV, rrie = 157 GeV; both models are marked in 
Fig. [2]) and compare it to the PAMELA data. For com- 
parison, we also show the expected background flux [2^ . 
Propagation effects thus considerably smear the spec- 
trum shown in Fig. [H but the clearly pronounced cutoff 
at Ef,+ = still remains as a prominent feature. It is 
interesting to note that this type of pronounced spectral 
signature so far has only been associated to Kaluza-Klein 
dark matter [sOl. Even though the cutoff in this latter 
case appears, due to the large branching ratio into e+e", 
to be even more pronounced, it would be observationally 
very challenging to see this difference with an energy res- 
olution of the about 5% expected for PAMELA. The ap- 
parent discrepancy between the background expectation 
and the new data at small energies is most likely due to 
a change in the solar potential which has not been taken 
into account so far [3|; this effect, however, is expected to 
be negligible at positron energies above around 10 GeV 
and we will therefore not discuss it further here. 
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FIG. 3: The solid line is the expected flux ratio /{e^ + e~) 
as calculated following [2^. The data points are the combined 
HEAT [ijl and PAMELA data Furthermore, the expected 
flux ratio for our benchmark models is shown without (dot- 
ted lines) and after taking into account radiative corrections 
(dashed lines). See text for further details. 



Unfortunately, the need for large boost factors is 
generic for all models that show a high positron yield en- 
hancement in the way reported here (in Fig. [31 we used 
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boost factors of 3 x 10^); the reason being that models 
in the coannihilation region have very small annihilation 
rates to start with. While such large boost factors are 
difficult to achieve in standard scenarios [s^l, they are 
easily encountered if, e.g., the Milky Way hosts a typical 
population of intermediate mass black holes [s^; other 
possibilities include a non-thermal neutralino production 
in the early universe (see, e.g., [31), a non-standard pre- 
BBN expansion rate [sst . or a very nearby dark matter 
"clump" (which, however, is quite unlikely according to 
present models of structure formation). 

One should also keep in mind that boost factors of at 
least 50-100 are needed in most cases, anyway [11], to see 
the effect of supersymmetric dark matter annihilation in 
the positron spectrum - but even for very large boosts, 
the resulting positron spectra are too soft to explain the 
observed steep rise in the jie^ ■\- e~) ratio. As be- 
comes apparent from Fig. [3l this is actually also true for 
the already quite hard spectra reported here in case of 
masses higher than rriy^ > 100 GeV. In order to really 
fit the PAMELA data through primary positrons from 
neutralino annihilations would thus require rather small 
neutralino masses. A generic prediction of this model is 
therefore that a sharp cutoff in the spectrum has to be 
observed already at energies only slightly higher than so 
far accessible. Such a well-pronounced, step-like feature 
would be a spectacular discovery in the next release of 
PAMELA data, or in future experiments like AMS-02. 

In concluding, let us briefly address the consequences 
of this type of radiative corrections for possible dark mat- 
ter induced contributions in other cosmic ray species. For 
neutrinos, for the same reasons as discussed for secondary 
positrons, the only chance for large effects appears in sit- 
uations with great enhancements of the annihilation into 
leptons, i.e. the channels /i"'"/x~7 and t~^t~j. Still, just 
as in the case of positrons, these channels are unlikely to 
have a large impact on the total neutrino yield. A po- 
tentially interesting source for antiprotons, on the other 
hand, are gluons from the annihilation into ttg final states 
- which, in the stop coannihilation region, is considerably 
enhanced compared to the lowest order result in exactly 
the same way as the ti^ channel discussed in [2l|. How- 
ever, since the lowest order annihilation rate is extremely 
small in the region of interest, the resulting absolute yield 
is still too small to be of great significance. This general 
expectation is in agreement with earlier studies [37| . 

To summarize, we have shown that radiative correc- 
tions may significantly enhance the dark-matter induced 
positron yield and result in a pronounced spectral signa- 
ture, a rising positron to electron ratio and a sharp cutoff 
in the positron spectrum at the neutralino mass m^. To 
obtain such a spectral feature, similar to that observed 
by PAMELA [i], very large boost factors are needed. On 
the other hand, if such a feature is observed, a strong en- 
hancement can also be expected in the gamma-ray flux 
at photon energies close to (2l| (while the impact on, 



e.g., the expected antiproton spectrum would be negli- 
gible); such a cross-correlation would of course provide 
even stronger evidence for the dark matter nature of the 
signal. An unambiguous, testable prediction of this class 
of models is that the positron excess will be cut off at 
an energy not too far from the maximal energy presently 
reported by the PAMELA collaboration, as larger masses 
do not reproduce the slope of the rising positron ratio (for 
a very similar situation for another class of dark matter 
particles, see [ss!]). 

Finally, let us mention that the radiative corrections 
to the positron yield from neutralino annihilations that 
have been reported here have been implemented in the 
current, publically available version 5.0.1 of DarkSUSY 
11. 
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